Elementary Excitations of Quantum Critical 2+1 D Antiferromagnets 
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It has been proposed that there are degrees of freedom intrinsic to quantum critical points that can 
contribute to quantum critical physics. We point out that this conclusion is quite general below the 
upper critical dimension. We show that in 2 + 1 D antiferromagnets skyrmion excitations are stable 
at criticality and identify them as the critical excitations. We found exact solutions composed of 
skyrmion and antiskyrmion superpositions, which we call topolons. We include the topolons in the 
partition function and renormalize by integrating out small size topolons and short wavelength spin 
waves. We obtain correlation length exponent v = 0.9297 and anomalous dimension rj — 0.3381. 
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Quantum phase transitions have been a subject of the- 
oretical and experimental exploration since the pioneer- 
ing work of John Hertz d]. Since then, quantum criti- 
cal behavior has been understood and studied as arising 
from quantum fluctuations of the order parameter [H [3] ■ 
In this traditional approach the quantum transition is 
studied via the Wilson rcnormalization group in which 
fluctuations of the order parameter are taken properly 
into account. It is said that quantum phase transitions 
follow the Landau- Ginzburg- Wilson paradigm (LGW). 

It has recently been suggested that quantum critical 
points will have properties that cannot be obtained from 
LGW order parameter fluctuations alone [3, 4]. In partic- 
ular, it was suggested that quantum critical points will 
have low energy elementary excitations intrinsic to the 
critical point whose fluctuations will contribute and mod- 
ify the critical properties. It was postulated that these 
excitations will be fractionalized[l, 

That some quantum critical points have elementary 
excitations different from those of each of the phases it 
separates can be inferred quite generally. We concentrate 
in relativistic quantum critical points, but we emphasize 
that these physics can take place in other systems. For 
such a system, which we take to be an antiferromagnet, 
we are interested in the Neel magnetization Green's func- 
tion, or staggered magnetic susceptibility. 

In the ordered phase the transverse Green function or 
susceptibility corresponds to spin wave propagation and 
it has a nonanalyticity in the form of a pole corresponding 
to such propagation: 
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Here Z{Lo,k^ is between and 1, and the incoherent back- 
ground Gincoh vanishes at long wavelengths and small 
frequencies. The fact that the Green's function has a 
pole means that transverse Goldstone spin waves are low 
energy eigenstates of the antiferromagnet. At criticality, 
the system has no Neel order and thus Goldstones cannot 
be elementary excitations of the system. 



In the disordered phase the Green function or suscep- 
tibility corresponds to spin wave propagation with all 
three polarizations and it has a pole nonanalyticity cor- 
responding to such propagation: 
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Here A{ijj, k) is between and 1, and the incoherent back- 
ground Gincoh vanishes at long wavelengths and small 
frequencies, A is the gap to excitations in the disordered 
phase. That this Green's function has a pole means that 
triplet or triplon spin waves are low energy eigenstates 
of the disordered antiferromagnet. For 2-1-1 D antiferro- 
magnets, and in general for antiferromagnets below the 
upper critical dimension, the quasiparticle pole residue 
A vanishes as the system is tuned to the quantum crit- 
ical point At criticality, triplon excitations have no 
spectral weight and thus triplons cannot be elementary 
excitations of the system. 

On the other hand right at criticality the response 
function below the upper critical dimension (below which 
77 7^ 0, while above 77 = 0) has nonanalyticities that are 
worse than poles 
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as obtained from the rcnormalization group studies of 
the nonlinear sigma model Q. Below the upper crit- 
ical dimension 77 is a nonintegral universal number for 
each dimensionality. This critical susceptibility has no 
pole structure, but has a branch cut. It sharply diverges 
at UJ = ck and is purely imaginary for lo > ck. Branch 
cuts in quantum many-body or field theory represent im- 
mediate decay of the quantity whose Green function is 
being evaluated. Hence the elementary excitations or 
eigenstates of the noncritical quantum mechanical phases 
break up as soon as they are produced when the system 
is tuned to criticality: they do not have integrity. The 
complete lack of pole structure and the branch cut singu- 
larity below the upper critical dimension mean that the 
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elementary excitations of the quantum mechanical phases 
away from criticality, the spin waves, cannot even be ap- 
proximate eigenstates at criticality as they are absolutely 
unstable. 

The quantum critical point is a unique quantum me- 
chanical phase of matter, which under any small per- 
turbation becomes one of the phases it separates. It is a 
repulsive fixed point of the renormalization group. As far 
as the transition from one quantum mechanical phase to 
the other is continuous, and both phases have different 
physical properties, the critical point will have its unique 
physical properties different from the phases it separates. 
The properties of the critical point follow from the critical 
Hamiltonian H{gc) {qc is the critical coupling constant), 
which will have a unique ground state and a collection 
of low energy eigenstates which are its elementary exci- 
tations. These low energy eigenstates are different from 
those of each of the phases as long as we are below the 
upper critical dimension. As a matter of principle, all 
quantum critical points below the upper critical dimen- 
sion will have their intrinsic elementary excitations. 

We have seen that below the upper critical dimension, 
the excitations of the stable quantum phases of the sys- 
tem become absolutely unstable and decay when the sys- 
tem is tuned to criticality. The question comes to mind 
immediately: what could they be decaying into? When 
one tries to create an elementary excitation of one of the 
phases, it will decay immediately into the elementary ex- 
citations of the critical point. The critical excitations will 
be bound states of the excitations of the stable phases 
the critical point separates. These bound states could be 
fractionalized as conjectured by Laughlin|3i] and Senthil, 
et. aZ. 4], but they need not be in all cases. These criti- 
cal degrees of freedom are responsible for corrections to 
the LGW phase transition canon [^]. The intrinsic quan- 
tum critical excitations contribute to the thermodynam- 
ical and/or physical properties of the quantum critical 
system. 

Now we turn to a specific model in order to identify 
what the critical excitations are. We concentrate on 2 -t- 
1 dimensional short-range Heisenberg antiferromagnets 
in a bipartite lattice. These are described by the 0(3) 
nonlinear sigma model augmented by Berry phases[7l@] 
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where ps = JS^ is the spin stiffness, and the spin-wave 
velocity c = 2y/2JSa with a the lattice constant and J 
the microscopic spin exchange. The Berry phase terms 
represent the sums of the areas swept by the vectors fti^r) 
on the surface of a unit sphere as they evolve in Euclidean 
time [7]. They were shown to be zero0, Q in the Neel 
phase and the critical point, but relevant in the disor- 



dered phase [lO[. We drop the Berry phase terms as we 
will study the critical properties of 2 -I- 1 D antiferromag- 
nets as approached from the Neel phase. 

The 0(3) nonlinear sigma model action ^ has a classi- 
cal "ground state" or lowest energy state with Neel order 
corresponding to a constant magnetization. The equa- 
tions of motion that follow from the action have approx- 
imate time dependent solutions, corresponding to Gold- 
stone spin wave excitations. The equations of motion, in 
2+1 D only, also have exact static solitonic solutions of 
finite energy [11,] 
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These solitons are called skvrmions[12j. 

Skyrmions are of a topological nature as they are char- 
acterized by the integer winding number 
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These configurations consist in the order parameter ro- 
tating an integer number of times as one moves from in- 
finity toward a fixed but arbitrary position in the plane. 
Since two dimensional space can be thought of as an in- 
finite 2 dimensional sphere where the magnetic moments 
live, the excitations fall in homotopy classes of a 2D 
sphere into a 2D sphere: S^]!^. Skyrmions ro- 

tate at finite length scales but relax into the Neel state 
far away: lim|5|^oo n = (0,0,-1). They have a direc- 
tionality given by the direction of the Neel order they 
relax to at infinity. The skyrmion number is a conserved 
quantum number as it is the zeroth component of the 
current = (l/87r)e^''°'n • 8,^71 x daii which can easily 
be checked to be conserved 5^ = 0. 

In order to study skyrmion properties more conve- 
niently, we use a very useful way of describing the 0(3) 
nonlinear sigma model, which is through the stereo- 
graphic proiectionflllfist: 
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In terms of w the nonlinear cr-model action is 
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where z = x + iy and z* = x ~ iy is its conjugate, and 
A/gA — Ps- gA is the microscopic Goldstone coiipling 
constant defined at the microscopic cutoff scale A[3, @]. 
The classical equations of motion which follow by sta- 
tionarity of the classical action are 



OqW — Adzdz'W 
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When the system Neel orders, n, or equivalently w, will 
acquire an expectation value: {n"') = S^"', (w) = 0. 

As mentioned above, besides Goldstones, there are 
static skyrmion configurations 11, w = Tli^i ^/{z — 
ai) whose topological invariant (jj]) in terms of the stere- 
ographic variable, w, is 
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This configuration can easily be checked to have charge 
q and energy AnqA/g^. y is the arbitrary size and 
phase of the configuration and ai are the positions of 
the skyrmions that constitute the multiskyrmion config- 
uration. Similarly, the multiantiskyrniion configuration 
can be shown to be w = HLi ^* l{z* — a*) with charge 
—q and energy ATrqA/g\. 

We now investigate whether skyrmions and anti- 
skyrmion configurations are relevant at the quantum 
critical point. As mentioned above, their classical 
energy is 4:TrA/g\, which is independent of the size of 
the skyrmion A. On the other hand, in real physical 
systems there are quantum and thermal fluctuations. 
These renormalize the effective coupling constant of the 
nonlinear sigma model and makes it scale dependent. 
To one loop order the rcnormalized coupling constant is 
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Since the skyrmion has an effective size A, spin waves of 
wavelength smaller than A renormalize the energy of the 
skyrmion via the coupling constant renormalization lead- 
ing to an energy and Euclidean action Se = I3E, which 
are now scale dependent through the scale dependence of 
the coupling constant at scale ji — 1/A. If the system is 
at temperature T = l/f3, this temperature sets the size of 
the skyrmion to be the thermal wavelength A = /3. The 
skyrmion Euclidean action is then 
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Having obtained the Euclidean action for skyrmions 
(H]), we now study its low temperature limit in the Neel 
ordered phase and at the quantum critical point. Ac- 
cording to the one loop renormalized coupling constant 
([5]), the quantum critical point occurs when the renor- 
malized spin stiffness (psifJ-) oc fJ-fg^,) vanishes at long 
wavelengths (/i — > 0): gA — gc = 27r^. Since the skyrmion 
gap is iTTjj,/ g^, the critical point corresponds to skyrmion 
gap collapse. When in the Neel ordered phase, gA < 9c, 
the skyrmion Euclidean action ([9]) is infinite. Therefore, 
the probability for skyrmion contributions is suppressed 
exponentially at low temperatures, vanishing at zero tem- 
perature. Skyrmions are gapped and hence irrelevant to 
low temperature physics in the Neel ordered phase. 



At the quantum critical point 5a 
skyrmion Euclidean action is 
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This action is finite and constant at all temperatures and 
in particular, it will have a nonzero limit as the tem- 
perature goes to 0: the skyrmion probability is nonzero 
and constant at arbitrarily low temperatures and zero 
temperature. Hence there are skyrmion excitations at 
criticality at arbitrarily low energies and temperatures, 
including zero at zero temperature. Therefore skyrmion 
excitations contribute to quantum critical physics. 

We have seen that skyrmions are relevant at critical- 
ity as the critical point is associated with skyrmion gap 
collapse and they have a nonzero probability to be ex- 
cited at arbitrarily low temperature at criticality. On 
the other hand, skyrmions have nonzero conserved topo- 
logical number while the ground state has zero skyrmion 
number. Absent any external sources that can couple di- 
rectly to skyrmion number, they will always be created in 
equal numbers of skyrmions and antiskyrmions. There- 
fore, in order to study the effect of skyrmions and anti- 
skyrmions we need to include configurations with equal 
number of skyrmions and antiskyrmions in the partition 
function or path integral. We found a time independent 
solution to the equations of motion given by 
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where A is the size of the configuration, 9 and ip are 
the arbitrary directions of the configurations, a is the ar- 
bitrary position of the configuration and n is an integer. 
This configuration is topologically trivial because it has 
q = as obtained from On the other hand, it is 
composed of arbitrary superpositions of equal numbers 
of skyrmions and antiskyrmions with q = ±n, i.e. the 
precise superpositions we need to sum over in the path 
integral for the q = sector. Since this configurations 
is made of topologically nontrivial skyrmions and anti- 
skyrmions, we dub it a topolon. While the argument of 
the tangent is obviously a sum of an n skyrmion and an n 
antiskyrmion, it appears to not be a fully general one as 
all the skyrmions are at the same position. By starting 
with a fully general skyrmion configuration and making a 
change of variables to an effective "center of mass" coor- 
dinate, it follows that the results are the same as having 
all skyrmions at the same place. 

The topolon with spatial and temporal size A has Eu- 
clidean action (161) 
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The partition function including topolon configurations 
is given by Z = J2n=o where 
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is the usual partition function for the nonhnear sigma 
model with no topolons and only the spin wave like fields 
v and Se is the Euclidean action for the nonlinear sigma 
model in terms of the stereographic projection variables 
We also have that 
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The Zn^Q is the path integral with the n topolons with 
spin waves v. Besides integrating over the spin wave con- 
figurations, we must integrate over the topolon parame- 
ters: its size A normalized to the lattice spacing 1/A, its 
position a normalized to the area A of the system, and 
over the solid angle of its orientation normalized to Att. 

In order to renormalize the nonlinear sigma model, 
we integrate the spin wave degrees of freedom i> with 
momenta between the microscopic cutoff A and a lower 
cutoff or renormalization scale /x[6|. We also integrate 
topolons of size between the microscopic minimum length 
1/A and a new larger renormalization length 1/fi. We 
thus obtain a double expansion in (1 — /i/A) and the cou- 
pling constant which leads to the renormalized action, 
the renormalized spin stiffness and the beta function 
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The last term in the spin stiffness and in the beta function 
is the contribution from the topolons and the rest is the 
contribution of the spin waves. The coupling constant at 
the quantum critical point gc — 23.0764 is obtained from 
Pi9c) = 0. 

The correlation length at scale /i is given by[6[ 
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The correlation length exponent is = —l/f3'{gc) = 
— (d(3 / dg\g^g^)^^ . The correlation length exponent with 
topolon contributions evaluates to v ~ 0.9297. The d = 
2 + e expansion of the 0{N) vector model, which agrees 
with the 1/A^ expansion for large TV, gives i/ = O.sfljj. 



We note that our value is larger than the accepted nu- 
merical evaluations of critical exponents in the Heisen- 
berg model, ly = 0.71125[3|, but about as close to this 
accepted Heisenberg value than the 2 + e expansion or 
the 1/N expansion. We conjecture that the difference 
between our value and the Heisenberg value is real and 
attributable to quantum critical degrees of freedom. 

Goldstone renormalizations of the ordering direction 
(7 = 713, and hence of the anomalous dimension 77, are no- 
toriously inaccurate. The one loop approximation leads 
to a value of = 2, thousands of percent different 
from the accepted numerical value of 77 ~ 0.0375[l5|. 
The large N approximation, which sums bubble dia- 
grams, is a lot more accurate. To order 1/7V one obtains 
rj = 8/(37r2iV) ~ 0.09 for iV = 3. We now calculate the 
value of rj from topological nontrivial configurations 
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For the anomalous dimension at the quantum critical 
point we obtain ij{gc) — 0.3381. On the other hand, we 
have seen that spin wave contributions tend to give quite 
large and nonsensical values of 77. In fact so large as to 
wash out the momentum dependence of the propagator. 
Hence, to calculate 77, spin wave contributions prove to 
be tough to control. Our calculation gives a value quite 
larger than the a ccep ted numerical value. We have re- 
cently calculated [17] the unique value of rj that follows 
from quantum critical fractionalization into spinous and 
find f] ~ 1. While our value obtained from topolons is 
far from 1, it is a lot closer than the accepted numerical 
Heisenberg value and the 1/N value. 
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